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We derive a microscopic effective action for superconducting contacts with arbitrary transmission 
distribution of conducting channels. Provided fluctuations of the Josephson phase remain sufficiently 
small our formalism allows to fully describe fluctuation and interaction effects in such systems. As 
compared to the well studied tunneling limit our analysis yields a number of qualitatively new fea- 
tures which occur due to the presence of subgap Andreev bound states in the system. We investigate 
the equilibrium supercurrent noise and evaluate the electron-electron interaction correction to the 
Josephson current across superconducting contacts. At T = this correction is found to vanish for 
fully transparent contacts indicating the absence of Coulomb effects in this limit. 



I. INTRODUCTION 

It is well known that supercurrent can flow through 
a non-superconducting barrier between two supercon- 
ducting reservoirs. Initially this effect was predicted^ 
and microscopically analyzed^ for a specific case of (usu- 
ally very thin) tunnel insulating barriers. Later it 
was understood that non-dissipative transport of Cooper 
pairs between two superconductors is also possible in 
many other types of weak links, such as, e.g., quan- 
tum point contacts^ and superconductor-normal-metal- 
superconductor (SNS) junctions^, i.e. if a piece of a 
normal metal is placed in-between two superconductors. 
In contrast to tunnel junctions, in SNS systems at suf- 
ficiently low temperatures appreciable supercurrent can 
flow even though a normal layer can be as thick as few 
microns. 

It turned out that the Josephson effect in supercon- 
ducting weak links without tunnel barriers is directly re- 
lated to another fundamentally important phenomenon: 
Andreev reflection^. Suffering Andreev reflections at 
both NS interfaces, quasiparticles with energies below 
the superconducting gap are effectively "trapped" inside 
the junction forming a discrete set of levels which can be 
tuned by passing the supercurrent across the system. At 
the same time, these subgap Andreev levels themselves 
contribute to the supercurrent thus making the behavior 
of superconducting point contacts and SNS junctions 
in many respects different from that of tunnel barriers. 
For an extended review summarizing various features of 
dc Josephson effect in different types of superconducting 
weak links we refer the reader to Refs. [ZHSl- 

The number of Cooper pairs transferred between two 
superconductors and, hence, the Josephson current can 
fluctuate around its mean valued—. While at non-zero 
temperatures thermal fluctuations of the supercurrent 
should naturally exist in all types of weak links, in the 
limit T — > the relevant physics is essentially determined 
by the presence or absence of subgap Andreev states. 
Provided such states are present fluctuations of the su- 
percurrent do in general occur even at T = and at 
subgap frequencies. E.g., the equilibrium supercurrent 



correlation functions show pronounced peaks at frequen- 
cies equal to the distance between Andreev levels inside 
the weak link. The amplitudes of such peaks turn out to 
scale asi£ J2 n ^n(l — ^n)j where T n is the normal trans- 
mission of the n-th conducting mode of the barrier and 
the sum is taken over all such modes. The latter de- 
pendence implies that ground state fluctuations of the 
supercurrent can be expected neither in the limit of low 
barrier transmissions T n — > (i.e. in tunnel barriers 
where no Andreev states arc present) nor in fully open 
contacts with T n 1. 



Note that the above considerations remain applicable 
if one can neglect Coulomb effects. In small-size super- 
conducting contacts, however, such effects can be impor- 
tant and should in general be taken into account. A lot is 
known about interplay between fluctuations and charging 
effects in superconducting tunnel barriers^. Here we ex- 
amine the properties of superconducting junctions going 
beyond the tunneling limit. We will analyze fluctuation 
and interaction effects and demonstrate that Coulomb 
blockade in such junctions weakens with increasing bar- 
rier transmissions and eventually disappears in the limit 
of fully open superconducting contacts. 



The structure of our paper is as follows. In Sec. II 
we derive an effective action for superconducting con- 
tacts with arbitrary distribution of channel transmissions 
which enables one to describe equilibrium fluctuations of 
the current and interaction effects. In Sec. Ill we make 
use of this action and evaluate the supercurrent noise 
in superconducting contacts. Low frequency current re- 
sponse and capacitance renormalization due to retarda- 
tion effects are discussed in Sec. IV. In Sec. V we ana- 
lyze the interaction correction to the Josephson current. 
A brief summary of our main results is presented in Sec. 
VI. Some general expressions and technical details are 
relegated to Appendix. 
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FIG. 1: Short coherent conductor between two superconduct- 
ing reservoirs. 



II. EFFECTIVE ACTION AND PHASE 
FLUCTUATIONS 



In what follows we will adopt the standard model of a 
superconducting contact and consider two big supercon- 
ductors connected with each other via a normal conduc- 
tor (see Fig. 1) characterized by arbitrary transmission 
distribution T n of its spin-degenerate conducting chan- 
nels. Below we will only consider the limit of sufficiently 
short normal conductors with effective Thouless energy 
£Th strongly exceeding the superconducting gap A in 
both reservoirs, exh 3> A. In addition, the normal con- 
ductor length is assumed to be much shorter than de- 
phasing and inelastic relaxation lengths. Coulomb inter- 
action between electrons in the contact area is described 
in a standard manner by an effective capacitance C . 

We will assume that the contact is biased by external 
current / which does not exceed the critical one Ic and, 
hence, can flow through the contact without any dissi- 
pation, i.e. I = Is- In the absence of fluctuations this 
external current sets the value of the order parameter 
phase difference \ — XL — XR between two superconduc- 
tors. The corresponding implicit dependence of x on the 
supercurrent Is has the form 3 



I s(x) 



x tanh 



_ eA sinx y Tn 

2 n yJl-T n sm 2 ( X /2) 

AyJl-T n sm 2 ( X /2) 
2T ' 



(1) 



Here and below we set h = 1 and define the electron 
charge to be — e. 

In order to analyze fluctuation and interaction effects 
in such superconducting contacts we will allow for fluctu- 
ations of the superconducting phase difference around its 
average value x an d employ the effective action formal- 
ism combined with the scattering matrix technique. This 
approach was proven to be very successful in the case 
of normal conductors^— and NS hybrid structure d 18 ' 19 . 
Following the usual procedure we express the kernel J of 
the evolution operator on the Keldysh contour in terms 
of a path integral over the fermionic fields which can be 
integrated out after the standard Hubbard-Stratonovich 
decoupling of the interacting termi^. Then the kernel J 



J = J VtpiVif2 exp(iS c [tp] + iSt[<p]), 



(2) 



where the terms S c {tp] and St[<p] account respectively 
for charging effects and for the transfer of electrons 
and Cooper pairs between two superconducting reser- 
voirs. Both these terms represent the functionals of the 
fluctuating phase variables ^1,2 (^) defined on the for- 
ward and backward parts of the Keldysh contour and 
related to fluctuating voltages V± t 2 across the conductor 
as 1^1,2 (i) = eVi.2- With the aid of the Josephson re- 
lation one trivially identifies the superconducting phase 
difference on two branches of the Keldysh contour as 
X + 2pi,a(t). 

The charging term is taken in the standard formic 



S c [cp] 



C_ 

2e* 



dt'{tp\ 



C_ 

^2 



dt(p + (p 



(3) 



where we also introduced "classical" and "quantum" 
parts of the phase, respectively ip + = + <^2)/2 and 
<y9_ = ipi — ip 2- The structure of the term St[<p] is 
the same as in the normal case, one should only re- 
place normal propagators by 2 x 2 Green- Gorkov ma- 
trix function a 13 i 20 ' 21 . The corresponding result can be 
expressed in the formSi 



Si 



(4) 



where Ql,r are 4x4 Green-Keldysh matrices of the left 
and right superconducting electrodes. The product of 
these matrices implies time convolution and curly brack- 
ets denote anticommutation. 

Without loss of generality we can set the electric po- 
tential (and, hence, fluctuating phases) of the right su- 
perconducting terminal equal to zero. Then the Green- 
Keldysh matrix of this electrode can be written in a sim- 
ple form 



Qr = a 



9 R 9 K 
q A 



(5) 



where g R ' A are retarded and advanced 2x2 matrix func- 
tions 



~r,a 



£T3 



R,A 



A 2 



(6) 



and g = g F — Fg is the Keldysh matrix, where 
F(e) = tanh(e/2T) is the Fourier transform of F(t) — 
—iJ '/ sinh[7rTi] and fj are the Pauli matrices. For sim- 
plicity we choose the order parameter A of the right 
superconductor real and, hence, we can set A = iA-f^. 
In order to properly account for analytic properties 
of the functions ^ R ' A here we keep an infinitesimally 
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small imaginary part i8 which allows to define £ R ' A = 
±sgn e Ve 2 - A 2 for |e| > A and £ R > A = i\/A 2 - e 2 for 
|e| < A. 

The Green-Keldysh matrix Ql of the left supercon- 
ducting electrode reads 

Q L {v){t,t') = CM + (t)Cg(t,t')£M-(t')C, (7) 

where we defined the matrices 



and 



M± = 

exp[±i(| + y>i(t))73] 

exp[±i(f + i p 2 (t))f 3 ] 



Substituting the above expressions for Ql and Qr into 
Eq. Q we arrive at the action which fully describes 
transfer of electrons and Cooper pairs to all orders in T n . 
In the case of tunnel barriers the channel transmissions 
remain small and one can expand St in powers of T n . 
Keeping the lowest order terms ~ T n of this expansion 
one recovers the well-known Ambegaokar-Eckern-Schdn 
(AES) action 13 . Here, however, we will go beyond the 
tunneling limit and analyze fluctuation effects at arbi- 
trary transmission values T n . 

To this end we will proceed similarly to Ref. [lj^ and 
introduce the matrix 



X = 1 - T n /2 + (T„/4) {Q L ,Q R } I 



v±=o- 



(8) 



As the action St vanishes for <P—(t) = one has 
Trln Xq = 0. Making use of this property we can identi- 
cally transform the action Q to 



5 t = -^Trln[l+A ^oA'] 



(9) 



where 



X' = l + (T n /A)({Q L ,Q R }-2)-X . (10) 
With the aid of the above expressions for Ql,r we obtain 



X^ = 



S)*-e* - J 7p L tanh rr E± s ( e ± e «) 



1 _ I (t+iS) 2 

— iii- 



(11) 

where ^ 2 = e 2 — A 2 and the subgap Andreev level inside 
the contact with energies ±e n (x) are defined in a usual 

way as 



e n ( X ) - A^/l-T n sin 2 ( X /2). 



(12) 



Now let us assume that fluctuating phases <p±(t) (or 
fluctuating voltages) at the junction are sufficiently small 



and perform regular expansion of the exact effective ac- 
tion in powers of these phases. Expanding Ql(<p) up to 
the second order in (p± (thus finding the matrix X') , from 
Eq. (0 we obtain 



iSt 



dt/IsMP-ffl + iSn-Si, (13) 



where the supercurrent I,s(x) is defined in Eq. ([T]) and 

t t 

S R = J dt' J dt"K{t' -t") V -(t')v+{t"), (14) 


t t 

Si = f dt' J dt"l(t' -t")ip-{t')ip-{t") (15) 


with both kernels lZ(t) and I{t) being real functions. The 
general expressions for these functions turn out to be 
somewhat lengthy and for this reason are presented in 
Appendix. Here we only emphasize some of the proper- 
ties of TZ(t) and I(t). 

To begin with, it is straightforward to verify that in the 
lowest order in barrier transmissions T n the result (|13[1 - 
(|T5|) reduces to the standard AES action^ for tunnel 
barriers in the limit of small phase fluctuations. Qualita- 
tively new features emerge in higher orders in T„ being 
directly related to the presence of subgap Andreev levels 
±e n (x) inside the contact. Consider, for instance, the 
kernel X(i) defined in Eq. (|ATT) . It can be split into three 
contributions of different physical origin 



l{t)=l 1 {t)+l 2 {t)+l 3 (t). 



(16) 



The first of these terms, T\(t), represents the subgap con- 
tribution due to discrete Andreev states. The Fourier 
transform of this term has the form (cf. the first line in 
Eq. 



T 2 sin 2 X 



2 e 2 ( X )cosh 2 (e„( X )/27y 

T 2 (l-r»)sin 4 ( X /2) 2 

-jj- [1 + tanh (e„(x)/2T)J 



x[5(oj-2e n ( X )) + S(oj + 2e n ( X ))[ 



(17) 



It is obvious that this contribution is not contained in the 
AES action at all. The general expression for the second 
term 1 2 (t) is defined by the second and third lines of Eq. 
(|Alj) . In the limit of small barrier transmissions this term 

scales as X2 oc Tn 2 and, hence, is not contained in the 
AES action either. This contribution can be interpreted 
as the " interference term" between subgap Andreev levels 
and quasiparticle states above the gap. In the low tem- 
perature limit T — > the Fourier transform of this term 
T-iu differs from zero only at sufficiently high frequencies 
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|w| > A + e n (x)- At higher temperatures T > e n (x); 
however, I 2 u> vanishes only for \oj\ < A — e n (x) an d re ~ 
mains non-zero otherwise. Finally, the third term X^{t) 
accounts for the contribution of quasiparticles with ener- 
gies above the gap. The Fourier transform of this term 
I^u, is defined by the fourth and fifth lines of Eq. (|Alj) . 
In the high frequency limit w> Aor for A — > this term 
reduces to the standard result for a normal conductor 



UJ UJ 
9 2B COttl ^ 



(18) 



where Rn is the normal contact resistance determined 
by the Landauer formula 



R N or ^ n 



(19) 



Turning now to the function 7Z(t) in Eq. (fl"4")) we note 
that its Fourier transform can be represented as 1Z U — 
1Z' U +i*R-u, where both 1Z' U and TZ'^ are real functions. The 
function TZ' U is even in uj while TZ'^ is an odd function of 
uj, thus implying that the function 7Z(t) is real. 

The functions lZ(t) and I(t) are not independent. For 
instance, the Fourier transform 1Z'^ is related to I u by 
means of the fluctuation-dissipation relation 



TZ" = 21 u tanh 



UJ 

2T' 



(20) 



The two functions IZ'^ and TZ" are in turn linked to each 
other by the causality principle: the function 1Z(t) should 
vanish for t < 0. The general expression for TZ' U (IA2|) and 
further details are presented in Appendix. 

Finally we would like to point out that with the aid of 
the above Gaussian effective action one can easily eval- 
uate the phase-phase correlation functions for our prob- 
lem. Combining Eqs. (TT3"l) - (|15l) with §S§ one finds 



(ip + (t 1 )(p + (t 2 )} 



(21) 



^l m ( __1 ) coth—e-^-^, 

2tt \Cuj 2 /e 2 + TZ uj J 2T 



(VJ+(ti)v-(t 2 )) = 

oo 

doj ( 1 



-iu)(ti —t 2 ) 



2ir \Coj 2 /e 2 +1Z u 

— OO 



Note that these expressions do not include the effect of 
(possibly existing) external impedance which we do not 
specify here. If needed, corresponding modifications can 
easily be implemented in a standard manner—. 

In addition to the above phase-phase correlation func- 
tions in what follows we will also need to define the ex- 
pectation value of the current operator 



(I(t)) 



V<p± 



JSc[<p]+iSt[<p] 



(22) 



and the current-current correlation function (I(t)I(t')). 
For the symmetrized version of this correlator we have 



(i(t)i(t') + f(t')i(t)) 



-2e< 



Sip-(t)8(p-(t') 



(23) 



III. EQUILIBRIUM SUPERCURRENT NOISE 

We first employ our results in order to describe fluctu- 
ation effects in superconducting contacts in the absence 
of electron-electron interactions. In this case after per- 
forming functional derivatives with respect to tp— (t) in 
Eqs. (|2"2"|) and (f2"3"| one should formally set ip_ = 0. 

Let us define the noise spectrum as 



1 



(J(ti)J(* 2 ) + 7(t 2 )J(ti))-ii = 



doj 
2tt' 



(24) 



Then from Eq. (J23J) for oj < 1/R N C we easily find 

5 W = 2e 2 2 LJ . (25) 

Together with Eq. (jAll) this result provides the complete 
expression for the equilibrium noise power spectrum in 
superconducting contacts with arbitrary distribution of 
channel transmissions T n . 

In the low temperature limit T — > and at subgap 
frequencies < uj < 2A the above general result reduces 
to the following expression 



Sui — 



_rp3/2 



7rA 2 T 2 (l-T n ) S in 4 ( X /2) 
l-T„sin 2 ( X /2) 

x i A(we„(x) - A 2 (l + cos 



S(u-2e n ( X )) 



X)) 



s (x)((^-e„(x)) 5 



i(x)) 



x^-£n(x)) 2 -A 2 % - A - 6„( X )) 



(26) 



where 8(t) is the Heaviside step function. Eq. (|2^|) 
demonstrates that the contribution of each transmission 
channel to the noise spectrum at has a narrow peak at 
uj = 2e n (x) while at higher frequencies uj > A + e n (x) 
continuous noise spectrum sets in. For even higher fre- 
quencies uj > 2A also quasiparticles with energies above 
the gap contribute to the noise spectrum and in the high 
frequency limit uj ^> A Eqs. (|2"5j) , (|A1|) reduce to the 
standard Nyquist expression for normal conductors 



UJ UJ 

coth — — . 

R N 2T 



(27) 



We also note that the expression presented in the first 
line of our Eq. (126|) matches with the result previously 
derived in Ref. lift 

Let us discuss some properties of the quantum low fre- 
quency current noise (|26[) in more details. We observe 
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0.5 1.0 1.5 2.0 2.5 3.0 

to/A 



(fT71) . This additional thermal noise peak was previously 
discussed in Refs. [TollTTI . 

Finally, we would like to point out that very recently 
a general analysis of persistent current noise in normal 
rings was developed 2 ^. Similarly to our present findings, 
this analysis demonstrates that persistent current noise 
spectrum has the form of sharp peaks occuring at zero 
frequency and at frequencies determined by the interlevel 
distances for quantum states with nonzero transition ma- 
trix elements. In the low temperature limit the zero fre- 
quency peak disappears while the peaks at non-zero ui 
persist down to T = 0. Essentially the same situation is 
observed in superconducting contacts analyzed here. 



FIG. 2: Low temperature noise spectrum for diffusive super- 
conducting contacts at x — 2.5, 1.5, 0.5 (from top to bottom). 



that essentially depends both on the channel trans- 
mission values T n and on the phase difference x- The 
amplitude of the peak at the frequency w = 2e n (x) in- 
creases with T n at small transmissions and decreases at 
higher T n vanishing in the limit of perfect channel trans- 
mission T n — > 1 except for a special point x = 7r in which 
case the contribution of a fully open channel reduces to 
the universal peak at zero frequency. Combining this 
peak with the continuous spectrum contribution, for a 
fully open single channel at T = and \ — n we obtain 



=Tre"A 2 6(Lo)+e 2 A\ 1 



A 2 



{uj-A). (28) 



In the case of many conducting channels with differ- 
ent T n narrow peaks originating from different channels 
occur at different frequencies and a smoother noise spec- 
trum is observed. An important example is a diffu- 
sive conductor characterized by the the so-called bimodal 
transmission distribution 



P(T n ) cx 



1 



(29) 



Averaging the result (|26|) with this transmission distribu- 
tion we arrive at the equilibrium zero temperature noise 
spectrum of diffusive superconducting contacts. The cor- 
responding results are displayed in Fig. 2 for different 
values of the phase difference x- The noise spectrum 
is zero for cj < 2A| cos(x/2)|, it increases with to at 
lu > 2A| cos(x/2)| reaching the maximum at 



A 



yjl + cosx+ v/(l + cosx)(17 + cosx) 



V2 



(30) 



and showing cusps at uj = A(l + | cos(x/2)|) and ui = 2A. 
In the limit cj ^> 2A the Nyquist noise (|2T|) is recovered. 

At non-zero temperatures there appear additional con- 
tributions to the noise spectrum. In particular, an extra 
peak at zero frequency emerges with the amplitude which 
depends on temperature as cx cosh _2 (e„(x)/2T), cf. Eq. 



IV. CAPACITANCE RENORMALIZATION 

Let us now take into account small voltage fluctuations 
V(t) across the contact. Since our present consideration 
is restricted to small fluctuations of the phase f+(t) = 
e J Q dt'V(t') <C 1, the constant in time part of the voltage 
should be equal to zero and the Fourier amplitude of its 
fluctuating part should obey the condition eV^ <C lu. 

Under these conditions the total current / across the 
superconducting contact takes the form 

I = I s (x) + -<P+(t)-e J dt'll(t-t')<p + {t')+5I(t), (31) 

where Is(x) is the supercurrent (fT]), the term involving 
C represents the displacement current, the 7?.-dependent 
term accounts for the retarded current response on the 
fluctuating voltage V(t) and SI(t) is the stochastic con- 
tribution to the current with the correlator (SISI)^ = S u 
studied in the previous section. Eq. (|3"Tj) represents the 
quasiclassical Langevin equation describing small fluctu- 
ations of the Josephson phase in superconducting con- 
tacts. 

Let us analyze the Fourier amplitude I u of the current 
in the limit of small temperatures and frequencies 



T,w<2 £rl (x). 



(32) 



We remark that the condition (f3"2l may yield paramet- 
rically different restrictions for weakly and highly trans- 
parent channels, cf. Eq. (fl"2"l) . In the limit ([32]) the noise 
term SI in Eq. f|3 1 j) vanishes, while the kernel 7Z can be 
expanded in lu up to ~ to 1 terms. Then we obtain 



I u = eAip^ 



^T»(cosx + T„sin 4 (x/2)) 

2^ 



n 



(l-T„sin 2 (x/2)) 



3/2 



C*(x) 2 

W tf +u . 



(33) 



The first term in this expression accounts for the shift 
of x in Eq. (fT|) by 2ip +ul - The renormalized capacitance 
C* (x) involved in the second term of Eq. (|33|) is defined 
as 



C*(x) =C + SC( X ), 



(34) 
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where in the limit T — > we have 

1 '2-(2-T n )sin 2 (x/2) 



(l-r„sin 2 ( X /2)) 



T„sin 4 ( X /2) 
1 2T„(T„-2)sin 2 ( X /2) 



(35) 



+5 + T„ 



2-2(l + 2r rt )sin 2 ( X /2) 
T„sin 4 ( X /2) 



Let us analyze some important limiting cases of the 
above general expression for 5C{x)- I n the tunneling 
limit T n <C 1 this result reduces to 



SC( X ) 



3tt 



32Ai? 



N 



1 



cosx 



(36) 



The first - x-independent - term describes the well 
known capacitance renormalization in Josephson tunnel 
junctions due to quasiparticle tunneling^. The second 
- x _ dependent - term (which originates from the so- 
called /3-terms in the AES action 1 ^) is usually neglected 
in the literature. This approximation is justified provided 
Coulomb effects are pronounced, phase fluctuations are 
strong and, hence, (cosx)x ~~ * 0- Here, however, we are 
dealing with small phase fluctuations in which case the 
X-dependent terms need to be fully accounted for. 

In the case of small Josephson phases x "C 7r Eq. ([33)) 
reduces to the universal expression 



SC 



16AR 



N 



(37) 



which remains applicable for any distribution of channel 
transmissions. Provided all channels are transparent, i.e. 
T n ~ 1, Eq. 435]) yields 



5C 



16AR N cos 4 (x/4) 



(38) 



in the region < x < for X n °t too close to tt. Starting 
from 7r — x ~ (1 — Tn) 1 / 5 the function 5C(x) deviates 
from Eq. (|38|) and tends to 



SC 



— Y- 

4A ^ (1 



1 



(1-T„)3/ 2 - 



(39) 



for x = tt- In the important case of diffusive contacts 
averaging of Eq. (|35[) with the bimodal transmission dis- 
tribution yields 



SC: 



1.05 



AR n (tt~x) 2 



(40) 



for 7r — x *C tt, while for small values of x we again re- 
produce Eq. ([57]). We observe, that the renormalized 
capacitance (|40|) for diffusive superconducting contacts 
diverges as the phase difference x approaches n. This be- 
havior is quite natural since (i) the contribution of almost 
fully open channels (|39[) becomes large in this limit and 



u 





T„=0.9 / 


/Diffusive - 












r„=o.i : 









0.5 1.0 1.5 2.0 2.5 3.0 

X 

FIG. 3: The capacitance 5C(x) normalized by Co = 
n /IGARm ■ As indicated in the plot, three curves corre- 
spond to uniform barriers with channel transmissions T n = 
0.1,0.4,0.9, while the fourth curve corresponds to diffusive 
contacts. 



(ii) many such channels are available in diffusive barri- 
ers. The behavior of the renormalized capacitance is also 
illustrated in Fig. 3. 

Eq. (f3"3")) also allows to determine the low temperature 
Josephson plasma frequency ujj of oscillations near the 
bottom of the Josephson potential well. We obtain 



e 2 A ^ T»cosx 

c*(x) V 



T 2 sin 4 (x/2) 



(l-T„sin 2 (x/2)) 3 / 2 



(41) 



where C*(x) is defined in Eqs. (|34[). (|35[). Strictly speak- 
ing, this expression applies only provided the condition 
(|3"2")l is fulfilled. However, qualitatively it remains valid 
also at ujj ~ 2e„(x) up to a prefactor of order one. With 
this in mind, below we will employ the above expressions 
also in this case. 

In the limit of large geometric capacitance of the junc- 
tion C 3> SC the capacitance renormalization can be 
neglected. In this case we have oj,j oc \J\fC, In many 
cases, however, geometric capacitance turns out to be 
negligibly small so that C* ~ SC. In such cases the 
combination loj /A depends only on x and on the barrier 
transmissions. This situation will be considered below. 

For small x *C it and T the Josephson current (fTJ) 
reduces to Is(x) = 7r ^x/(2ei?Ar) for any transmission 
distribution. Combining this expression with Eq. (|37[) 
we get ujj/A — 4. This result universally holds for small 
values of the Josephson phase. For higher values of x the 
Josephson plasma frequency becomes smaller. E.g. in 
the case of tunnel barriers T n -C 1 for — 7r/2 < x < 7r /2 
one trivially finds 



¥± = 4 / 2cos X 
A V 3 — cos x 



(42) 



In the case of the highly transparent contacts the 
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FIG. 4: The Josephson plasma frequencies of highly trans- 
parent (T n = 0.9), diffusive and tunnel barriers (from top to 
bottom) in the limit SC 3> C. 



Josephson plasma frequency can be written as 



w / o X 

-J-=4cos 2 4 
A 4 



cosf, 



(43) 



where we assume that — ir < x < tt and 7r — |x| > (1 — 
Tn) 1 / 5 . As in this case the critical current is achieved at 
Xmax ~ 7r - 2(1 - T^) 1 / 4 , for x close to x-max we obtain 



^=2V2(1 



m«i X- 



(44) 



Provided the current is close to the critical one, the 
Josephson plasma frequency tends to zero as 



A 



7 



1 



J 2 

T'2 
1 C 



1/4 



(45) 



where 7 = 4^/2/3 for tunnel barriers, 7 ~ 2.53 in the 
diffusive limit and 7 = 2\/2(i — Tii) 1 ^ 8 f° r highly trans- 
parent junctions. The behavior of ojj{x) for tunnel, dif- 
fusive and highly transparent barriers is also illustrated 
in Fig. 4. 

Finally, let us display the adiabatic form of effective 
action for a superconducting contact which remains ap- 
plicable for small phase fluctuations under the condition 
(1321). It reads 



(46) 



dt' 



c*{x) 



-<P <P 



I s (x + ^+(t')) 



where Is(x) 1S defined in Eq. (J}. 



V. INTERACTION CORRECTION TO 
SUPERCURRENT 

Let us now turn to the electron-electron interaction 
correction to the equilibrium Josephson current ([T]). Pre- 
viously such correction was analyzed in the case of 



Josephson tunnel barriers in the presence of linear Ohmic 
dissipation^ (see alsoi^). The task at hand is to investi- 
gate the interaction correction to the supercurrent in con- 
tacts with arbitrary transmission distribution. As before 
in this paper, we do not include any external impedance 
into our consideration. 

In order to evaluate the interaction correction it is nec- 
essary to go beyond the Gaussian effective action (fTc 



(|TlM) and to evaluate the higher order contribution ~ <p 6 . 
It is easy to observe that the interaction correction to the 
supercurrent is provided by the following non-Gaussian 
terms in the effective action: 



S(iS t ) = 

J J J dhdt 2 dt 3 Y{ti 



(47) 



dt 1 dt 2 dt 3 Z(t ll t 2 ,t 3 )tp+ {ti)ip- (t 2 )cp- (t 3 ). 



The function Y(ti,t 2 ,t 2 ) can be written as 



Y(h,t 2 ,t 3 ) = 
dcui dcu 2 



(48) 



y(w 1 ,w 2 ) e - l " i ( ti ~ t2 ' e -^ 2( * i -* 3) , 



where Y(loi,uj 2 ) — Y(u> 2 ,uji). The function Zitx^.t^ 
can be expressed in a similar way. 

Adding the non-Gaussian terms ([47|) to the action and 
employing Eq. (|22l) we arrive at the following expression 
for the interaction correction 



SIsix) 



+2ie 



did 

2^ 



■^■Y(u),-u}){(p+(p + ) ( , 
Zn 



Z(0,-u)((p + ip_) u , 



(49) 



where the phase-phase correlators are defined in Eq. (|21[) . 

Let us consider the first term in the right-hand side of 
Eq. (|49|) . It is easy to see that in the limit of low tem- 
peratures only frequencies \ui\ > A + e„(x) contribute to 
the integral in Eq. (l2~Tj) for (ip + (p+) while the contribu- 
tion from the frequency interval \ui\ < A + e n (x) vanishes. 
Furthermore, the leading contribution from the first term 
in Eq. (|49f is picked up logarithmically from the interval 
2A < \u>\ < l/R N C where 



e 2 R N 



(50) 



and the function Y(uj, —uj) tends to a frequency indepen- 
dent value. 

After a straightforward but tedious calculation (some 
relevant details are presented in Appendix) in the inter- 
esting frequency range cj A from Eq. @ one finds 

iA sin% 



Y(u, 



(51) 



r ra (l-r n )(2-T n sin 2 ( X /2)) 
(l-T n sin 2 ( X /2))3/2 
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Similarly to the third term in Eq. (IA7[) this high- 
frequency term involves the factor 1 — T n , i.e. it van- 
ishes for fully open conducting channels. Combining Eqs. 
(f5Uj) . (jSTj) with we arrive at the expression for the 
supercurrent 



I(x) = Is(x)+SIs(x)- 



(52) 



In the limit of low temperatures the interaction correc- 
tion reads 



2gN 

T„(l 



In 



SIs(x) = 



l 



2AR N C 



(l-T„sin 2 ( X /2))3/ 



;( 2 



smx 



-T„ sin - 



(53) 



where = 2ir / (e 2 Rn) is the dimensionless normal state 
conductance of the contact. This result is justified as long 
as the Coulomb correction 51 s (x) remains much smaller 
than the non- interacting term Is{x) ©■ Typically this 
condition requires the dimensionless conductance to be 
large g N > ln(l/2A J R A rC). 

Note that Eq. (I5B1 was derived only from the first term 
in Eq. (|49[) . The second term in this equation involving 
the function Z(0, — uj) and the correlator (cp+tp-) can be 
treated analogously. The corresponding analysis demon- 
strates that the contribution of this term turns out to 
be smaller than that of the first term by the logarithmic 
factor ~ ln(l/2Ai?jvC)- Accordingly, the second term in 
Eq. (pi!)|) can be safely neglected for our purposes. 

Let us emphasize again an important property of the 
result (j53"|) : The interaction correction contains the fac- 
tor 1 — T n and, hence, vanishes for fully open barriers. In 
other words, no Coulomb blockade of the Josephson cur- 
rent is expected in fully transparent superconducting con- 
tacts. Note that this conclusion is also consistent with 
numerical results in Ref. [2~3 . 

The expression for the interaction correction (|53[) can 
further be specified in the case of diffusive contacts. In 
the absence of interactions the Josephson current in such 
contacts follows from (JXJ) and takes the well known form 



I six) 



TTA 



X, 1 
cos — m — 

2eR N 2 1 - sin | 



sin I 



(54) 



Including interactions and averaging (|53|) with the bi- 
modal transmission distribution (1291) one finds 



6Is(x) 



-A In 



2AR N C 



cot( X /2) (55) 



^sii 



sin — In — 
2/ 1 



sin(x/2) 
sin(x/2) 



Note that the result (|5"3")) can formally be reproduced 
if one substitutes T n — > T n + ST n into Eq. (fTJ , where 



ST n = — — In ( — — t 

g N \2AR N C 



T„(l - T n ), (56) 



and then expands the result to the first order in 
5T n . Interestingly, the same transmission renormaliza- 
tion (|56|) follows from the renormalization group (RG) 
equations J&il 



dL 



T n (l-T n ) 



L = In 



1 



(57) 



derived for normal conductors. In order to arrive at Eq. 
(1551) one should just start the RG flow at e — 
and stop it at e = 2A. Thus, the result ([53]) can be in- 
terpreted in a very simple manner: Coulomb interaction 
provides high frequency renormalization T n + ST n l|56[) 
of the barrier transmissions which should be substituted 
into the classical expression for the supercurrent ([TJ. It 
should be stressed, however, that the last step would by 
no means appear obvious without our rigorous derivation 
since the Coulomb correction to the Josephson current 
originates from the term ~ f>~^>\ in the effective action 
which is, of course, totally absent in the normal case. 



VI. DISCUSSION 

The analysis employed in this paper demonstrates 
that fluctuation and interaction effects in superconduct- 
ing contacts with arbitrary transmissions of conducting 
modes show a number of qualitatively new features as 
compared to the case of Josephson tunnel barriers^. The 
main physical reason behind such differences is the pres- 
ence of subgap Andreev bound states (fl"2"]) inside the sys- 
tem. 

In the limit of sufficiently small fluctuations of the 
Josephson phase difference we derived the complete ex- 
pression for the effective action of superconducting con- 
tacts. This expression allowed us to obtain the general re- 
sult for the equilibrium current-current correlation func- 
tion describing Josephson current noise in such contacts. 
Due to the presence of subgap bound states this cur- 
rent noise essentially depends on the Josephson phase x 
and remains non-zero even in the zero temperature limit 
and at subgap frequencies. For instance, in a physically 
important case of diffusive contacts at T = the equi- 
librium current noise spectrum differs from zero at all 
frequencies exceeding the threshold value 2A|cos(x/2)| 
and has the form displayed in Fig. 2. 

Another important effect studied here is capacitance 
renormalization. While geometric capacitance of super- 
conducting contacts can be small, retardation effects 
yield additional "capacitance-like" contributions which 
depend on both the transmission distribution T n and the 
Josephson phase x and can well exceed the geometric ca- 
pacitance term. In particular, in the limit x — ^ 77 the 
renormalized capacitance was found to diverge in highly 
transparent and diffusive contacts, see Fig. 3. This be- 
havior differs from that of tunnel junctions^ and can also 
substantially affect the frequency of Josephson plasma 
oscillations, see Fig. 4. 
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Finally, our effective action formalism enabled us to 
analyze the correction to the Josephson current due to 
electron-electron interactions. Provided this Coulomb 
correction Sic remains small as compared to the Joseph- 
son critical current Ic one finds the universal result 

f^ = -^2Aib)' (58) 

where the numerical prefactor a depends on the trans- 
mission distribution. This prefactor reaches its maximum 
value a = 2 in the case of tunnel barriers and becomes 
smaller for higher transmissions, e.g. a ~ 0.72 for diffu- 
sive contacts. In the case of fully open barriers the pref- 
actor a tends to zero, a — > 0, implying that no Coulomb 
blockade of the Josephson current is expected in such 
barriers. 

Our results for the interaction correction, e.g., Eq. 
(|58[) . might explain a rapid change between super- 
conducting and insulating behavior recently observed 25 
in comparatively short metallic wires with resistances 
close to the quantum resistance unit ~ 6.5 Kf2 in- 
between two bulk superconductors. Previously it was 
already argued 2 ^ that such a superconductor-to-insulator 
crossover can be due to Coulomb effects. Our present re- 



sults provide further quantitative arguments in favor of 
this conclusion. 

Let us again stress that in this paper we only consid- 
ered the limit of relatively short contacts in which case 
the characteristic Thouless energy of the contact exh ex- 
ceeds the superconducting order parameter A. In the 
opposite case of long junctions exh *C A (which appears 
to be more relevant to the experiments^ 5 -) the interaction 
correction to the supercurrent turns out to have the same 
structure ([55)1 with 2 A substituted by exh (see also 1 ^). 
The corresponding analysis will be published elsewhere. 
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Appendix A 

Let us present our results for the kernels X(t) and lZ(t) 
in Eqs. |[T3|). ([15]) . The Fourier transform 1^ of the kernel 
I(t) is even in Li, real and non-negative function. Hence, 
l(t) is a real function. For u > we have 



X - = E { A ' sin2 X [1 - F 2 (e n ( X ))} 5(c) + ^|iL-pA 4 sin 4 | [l + F 2 ( e „( X ))] S ( w - 2e n ( X )) (Al) 

+T 3 / 2 [1 - F(e n ( X ))F(u + e n ( X ))] 9(cj - A + e„( X )) 



+T„ 3 / 2 [1 + F(e n (x))F(u - e n ( X ))} 6(u - A - e n { X )) 
" de Ve 2 -A 2 v /(e + ^) 2 -A 2 



sin ■ 



. X 
sin — 
2 



A(^ K ( X )+A 2 (l + COS X )) n - - 

4e„(x) ((w + e„(x)) 2 - e 2 (x)) 

A(^„(x)-A 2 (l + cosx)) n . 2 

4e„(x) ((w - e„(x)) 2 - e 2 (x)) 



27T ( e 2- e 2 W )(( e + a; )2_ e 2 (x)) 
A 



e{oj + e) + A 2 cos x + T„ A 2 sin 2 J ) (1 - F(e)F(e + w)) 



^_2A1 / de V^A^V(^-g) 2 "A 2 
2 1 W 2^( e 2 - e 2 (x))((^- e ) 2 -£ 2 (x)) 



e(w - e) - A 2 cos x - T n A 2 sin 2 J ) (1 + F(e)F(w - e)) } , 



The Fourier transform of the function TZ(t) can be writ- 
ten in the form 1Z U = 1Z' U + iTZ^, where the real functions 
1Z' U and IZ'u are respectively even and odd in u>. Hence, 



lZ(t) is also a real function. The function 1Z" is linked 
to I by means of FDT relation(|20"]) , and the function 1Z' U 
for lu > reads 
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F (zn{x)) A 2„ ; „2 X ^(l-T n )F(e„(x)) A4s . n4 x 



Sm ~2 2 "(^„ 4e 2 (x))en(x) 



(A2) 



7 2tt (e 2 -e 2 (x))((e + ^) 2 ™e 2 (x)) 

A-min{cj,2A} 

, T„ 3/2 A|sin(x/2)|F(6„(x)) 



2e«(x) 
+0(A + e„(x)-w) 



0(A - e n ( X ) - w) 



A 2 (l + cosx)-we n (x) 



A 2 (l + cosx)+^6n(x) 
(w + e„(x)) 2 -e 2 l (x) 



VA 2 -( W + e„(x)) 2 + 



( w - e „(x)) 2 - £ 2 (x) 



VA2-( w -e„(x)) 2 



The integral in this expression is taken in the v. p. sense. 
The functions TZ' U and TZ'£ are also related by the causal- 
ity principle: the function lZ(t) should vanish for t < 0. 
E.g., the contribution of Andreev bound states to 7Z'^ has 
the form 



<(1 - T„)^|iMA 4 sin 4 | [J( W - 2 £ „(x)) 



-S(Lu + 2e n ( X ))}. 
It contributes to lZ(t) as 



(A3) 



,^Mx)) ,4-4 X . 



4(x) 



A 4 sin 4 ^sin[2e„(x)t]. (A4) 



The corresponding contribution from discrete Andreev 
states to 1Z' reads 



T 2 (l-T n )F(e n ( X )) 
( W 2-4e2(x))e„(x) 



A 4 sin 



;„4 X 



(A5) 



One half of this contribution comes from the second term 
in Eq. (IA2[) . the other half stems from the second Heav- 



iside function in Eq. (|A2[) . Making use of the integral 

oc 

cos(bx) 7r 



dx- 



x 2 — y 2 2y 



= -^-sin(|6|y), 



(A6) 



taken in the v. p. sense, we observe, that the contribution 
(IA4[) increases by the factor two for t > and it vanishes 
for t < 0. Hence, the total contribution of Andreev levels 
to TZ(t) is given by the last term of Eq. (|A7|) . 
One can also write 

f 2 r°° ( T 

Tl(t) = 22 l"- 6 *^) J du3sm(wt) [f u + -^o 

-*»5\t) +T„(1 - T„)^%^A 2 sin 2 ^(t) (A7) 
+2T 2 (l-^) F(e " (x)) ^ 4 -" 4X 



e 2 (x) 



A 4 sin 4 ^(i)sin(2e n (x)t) , 



where 



oo 

= /-de Ve^Ay(e + i.) 2 -A 2 

^ 'V 2^( £ 2 - e 2 (x) ) ((e + w) 2_ e 2 (x)) 

A 

m— A 

-T„0(w > 2A) 



e(w + e) + A 2 cosx + T„A 2 sin 2 20 (F(e) - F(e + w)) (A8) 



de Ve 2 -AV(w-e) ; 



T 3/2 A 
2e„(x) 



. X 
sm — 
2 



27T ( e 2_ e 2 (x))(([J _ e) 2_ e 2 (x)) 

we„(x) + A 2 (l + cosx) 



(e(w - e) - A 2 cosx - T n A 2 sin 2 |) (F(e) + F(w - e)) 



L (^ + e«(x)) 2 -4(x) 



V(e n (x)+^) 2 -A 2 + e„( X )) - ^M*))] ^ - A + e„(x)) 



+ T (X) ; A w + ^ - e "(x)) 2 - A 2 [F(u - c(x)) + ^n(x))] f(« - A - e„(x)) 



The second term in Eq. (|A7|) describes Ohmic damping correction. Note that this correction turns out to be pro- 
and the third contribution represents the supercurrent 
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portional to the factor T„(l — T n ). Finally, the last term 
accounts for the effect of discrete Andreev levels as it is 
described above. 

Setting e m = max{A,T}, for e m t <C 1 with the loga- 
rithmic accuracy we obtain 



(A9) 



9(t) / dujsm(ujt) f L 



-6{t)T n (A 2 cos x + T n A 2 sin 2 In 
7r \ Z I te m 

Let us also present some details relevant for the calcu- 
lation of the kernel Y(lu, — 10) in the limit lu ^> A. This 
kernel can be represented as a sum 



Y = Y (r) +Y {2) +F (3) 



(A10) 



where F^' 2 ' 3 ) dehne respectively the first, second and 
third order terms in the expansion of the logarithm in 
Eq. dU). They read 



Y 



(i) = iT n F(e n (x)) 
2 e n ( X ) 



A sinx, 



(All) 



y(2) 



iTZF(e n ( X )) A2 



e«(x) 



A^ sin x 



iTl F(e n ( X )) 
8 eUx) 



A 2 sinx(e 2 t (x)- A 2 cosx), 



(3) ^_ ! 7|FMx)) 4 X 
8 eUx) X 2 



Combining these expressions with Eq. (IA10[) we arrive 
at Eq. (I5T1) . Consider, e.g., a typical summand 



5{iSt)= 256 / / / *i*2*3^-(ti)v'+(t2)^ + (t3)Tr[(OX - 1 Q ii ) tiit Jf3Q i X - 1 gfl) t2)t Jf3Q i X - 1 Q il Q i ) t3 ^ 



{OQ L Q R X^Q R ) ti t2 {f z Q L X^Q R ) t2M (f 3 Q L X { 



>t 3 M 



(A12) 



which emerges in the third order. Here the matrices Ql,r are taken for tp± — but for non-zero X; the indices stand 
for the temporal arguments of the convolutions and 



- I o a ° I " 



f 3 
r 3 



(A13) 



The contribution of Eq. (|A12j) to Y^ (u, -u) reads 
5YW(u),-cj) = |g / ^Tr [(QlAo-^aQlOAo- 1 ^, - Q^o"^^^^)^ Mi^fe) J ■ ( Al4 ) 



r 



Taking the limit u) 3> A and using the property 

S ( X " ^ ( l^iSf^l? + (x + i<5) 2 - b 2 
= -^(«-6). 



(A15) 



we find 



5F(3) = A2sinx(A2c ° sx " e " (x)) 



(A16) 



Depending on the way of counting there are about 50 
contributions to iSt similar to Eq. (|A12I) . The calcula- 
tion of the corresponding contributions to Y is analogous 
to that presented above. 
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